Bicocca-FT-02-14 July 2002 



Critical parameters and universal amplitude ratios 
of two-dimensional spin-^S Ising models 
using high- and low-temperature expansions 

P. Butera* and M. Comi 

Istituto Nazionale di Fisica Nucleare and 
Dipartimento di Fisica, Universitd di Milano-Bicocca 
3 Piazza della Scienza, 20126 Milano, Italy 

A. J. Guttmannt 

Department of Mathematics and Statistics, 
The University of Melbourne, Victoria, 3010, Australia. 
(Dated: February 7, 2008) 

Abstract 

For the study of Ising models of general spin S on the square lattice, we have combined our 
recently extended high-temperature expansions with the low-temperature expansions derived some 
time ago by Enting, Guttmann and Jensen. We have computed for the first time various critical 
parameters and improved the estimates of others. Moreover the properties of hyperscaling and 
of universality (spin S independence) of exponents and of various dimensionless amplitude com- 
binations have been verified accurately. Assuming the validity of the lattice-lattice scaling, from 
our estimates of critical amplitudes for the square lattice we have also obtained estimates of the 
corresponding amplitudes for the spin S Ising model on the triangular, honeycomb, and kagome 
lattices. 
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I. INTRODUCTION 



The properties of the spin S — 1/2 two-dimensional Ising model with nearest-neighbor 
interactions in zero magnetic field, have been extensively explored in the last six decades. 
Much more modest efforts have been devoted to the study of the simplest generalizations of 
the model to spin S > 1/2. The main reason is probably that these models are not known 
to be solvable or, at least, to have any simple duality property which can help to extend 
the small body of information coming from numerical methods of limited accuracy such as 
stochastic simulations, series expansions or transfer-matrix calculations. 

The first important result from a comparative study of Ising models for different val- 
ues of the spin came from pioneering work by Domb and Sykes 1 . They analysed the 
high-temperature(HT) expansion of the susceptibility x{K\ through 0(K 6 ) in the three- 
dimensional case and conjectured that the value of the critical exponent j(S) is independent 
of the spin magnitude. This was the first step towards the modern formulation of the critical- 
universality hypothesis. Similar analyses were soon repeated by other authors using both 
HT 2 ' 3 and low-temperature (LT) expansions 4 for two-dimensional systems. Unfortunately, 
the series derived in those years were rather short and therefore the results of the anal- 
yses could not reach a sufficient accuracy or were inconclusive. It was only in 1980 that 
Nickel 5 ' 6 finally extended through 0(K 21 ) the HT series in two dimensions on the square 
(sq) lattice and in three dimensions on the body-centered cubic lattice. The expansions of 
x(K; S) and of the second moment of the spin-spin correlation function ^{K; S) were then 
published only for S — 1/2, 1, 2, oo. More recently also the LT expansions on the sq lattice 
for S = 1/2, 1, 3/2, 2, 5/2, 3 were considerably extended by Enting, Jensen and Guttmann 7 . 
We have summarized in Table I and Table II the state of HT 3 ' 8 and LT 7 expansions before 
our work. 

In spite of the very large number of LT expansion coefficients now available, the analysis 
of the series remains arduous due to occurrence 7 of numerous unphysical singularities in the 
complex temperature plane 9 which are closer to the origin than the physical singularity and 
whose structure becomes increasingly complicated with S. As a consequence, the LT study 
of Ref. 7 has been an alarming lesson on the subtleties in the analysis of slowly convergent 
series more than a source of accurate estimates of the critical parameters of the models. 

Many intriguing indications and conjectures about the structure of these unphysical sin- 
gularities also came in the same period from work by Matveev and Shrock 10 who examined 

TABLE I: The longest HT expansions, published (or obtainable from data in the literature) before 
our work 11 ' 12 , for the susceptibility x(if; S), the second moment of the correlation function H2(K\ S) 
and the second field-derivative of the susceptiblity X4(if; S) in the case of the Ising models with 
general spin S on two-dimensional lattices. It should be noted that in the special case S = 1/2, on 
the sq lattice, much longer expansions 13 for x and \i 2 have been computed. However, the published 
expansions 5 of Xi{K; 1/2) on the sq lattice do not extend beyond K 17 . 



Observable 


Lattice 


Order 


Ref. 


X (K;S), fi 2 (K;S) 


sq 


21 


6 


Xi{K-S) 


sq 


10 


8 


X (K;S), ,2 2 (K;S) 


tr 


10 


3 


Xi{K-S) 


tr 


10 


8 
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TABLE II: The longest LT expansions, presently available for the spontaneous magnetization, the 
specific heat and the susceptibility in the case of the Ising models with various values of the spin 



on the sq lattice. 

Observable Lattice Order Ref. 

M(u;l), C(u;l), x(u; 1) ^ 113 ~ 

Af(u;3/2), C(u;3/2), x(u;3/2) sq 100 7 

M(u;2), C(u;2), X (u;2) sq 119 7 

M(«;5/2), C(«;5/2), X («;5/2) sq 126 7 

M(u;3), C(u;3), x(«;3) sq 154 7 



the spin S models on various two-dimensional lattices using transfer-matrix methods. 

Here we discuss some results of an analysis of HT series for the sq lattice recently 
extended 11 by linked-cluster expansion techniques. 

For the nearest-neighbor correlation function G(K;S), for x{K',S), and fj,2(K;S) our 
series reach order K 25 , while for the second field derivative of the susceptibility Xti-K'i they 
extend through 0(K 23 ). In order to make alternative analyses possible, our vast collection of 
series data both for two- and three-dimensional lattices was made easily accessible 12 on-line 
for S = 1/2, 1, 3/2, 2, 5/2, 3, 7/2, 4, 5, oo. It should be noted that HT and LT expansions as 
extensive as those obtained by Nickel et al. in Ref. 5 and more recently in Ref. 13 (only for 
x) in the very special case of the (partially solvable) two-dimensional S — 1/2 Ising model 
seem presently beyond reach for S > 1/2. 

The HT series show somewhat simpler and faster convergence properties than the LT 
series, because the behavior of the coefficients is dominated by the physical singularity. 
Although, even in this case, these favorable properties slightly deteriorate for S > 2, we can 
hope to determine basic HT critical parameters with a reasonable accuracy for various values 
of S. Moreover it is also worthwhile to reconsider the LT expansions of Ref. 7 for the sq lattice, 
because by relying on the results of our HT analysis, we can improve some estimates of the 
LT critical parameters and thus obtain new determinations of universal combinations 14 of 
LT and HT amplitudes. No theoretical surprises are expected from this analysis, however we 
believe it is still useful to improve the rather modest numerical precision presently available 
even for basic critical parameters like the critical temperatures, to determine various critical 
amplitudes for which no estimates are yet known and to use our results to test with higher 
accuracy the validity of hyperscaling and of universality with respect to the magnitude of 
the spin. 

Almost all the computational effort in extending series for the two-dimensional Ising 
model for S > 1/2 has been devoted to square-lattice series. However by making use of the 
theory of lattice-lattice scaling, as developed by Betts et al. 15 and extended by Gaunt and 
Guttmann 16 , using our estimates of the critical amplitudes on the square lattice, we are able 
to calculate the corresponding amplitudes on other two-dimensional lattices to precisely the 
same precision as they are known for the square lattice. 

II. THE SPIN-S ISING MODELS 

The spin-S* Ising models with nearest-neighbor interaction are defined by the Hamiltonian: 
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H{s} = - J - £ s{x)s^)-h^s{x) (1) 

(a;, a;') x 

where J is the exchange coupling, and s(x) = s z (x)/S with s z (x) a classical spin variable 
at the lattice site x, taking the 2S + 1 values —S, —S + l, . . . , S — 1, S. The sum runs over all 
nearest-neighbor pairs of sites. We shall restrict ourselves to the square lattice and consider 
expansions either in the usual HT variable K = J/ksT and in the natural LT variable 
u(S) = exp(—K/S 2 ). Here T is the temperature, ks the Boltzmann constant, and K will 
be called "inverse temperature" for brevity. In the critical region we shall also refer to the 
standard reduced-temperature variable t(S) = T/T C (S) — 1 = K C (S)/K — 1. 

In the HT phase, the basic observables are the connected 2n-spin correlation functions. 
Our series 12 cover quantities related to the two-spin correlation functions (s(x)s(y)) c and to 
the four-spin correlation functions {s(x)s(y)s(z)s(i)) c . 

In the LT phase the symmetry is broken and the n-spin correlations are non trivial also 
for odd n. In particular, we shall reconsider the LT expansions of the magnetization, the 
susceptibility and the specific heat derived for S = 1, 3/2, 2, 5/2, 3 in Ref. 7 . 

The spontaneous magnetization is defined by 

M(T;S)= Inn <s(0)>. (2) 
The internal energy per spin is given in terms of the nearest-neighbor correlation function 

by 

U(T;S) = -f(s(0)s(5)) (3) 

where 5 is a nearest-neighbor lattice vector and q is the lattice coordination number. 

The specific heat is the temperature-derivative of the internal energy at fixed zero external 
field 

CH( r ;5 , = ^p. (4) 

In terms of x(T; S), the zero- field reduced susceptibility, 

X (T;<?) = X>(0) S (£)) c (5) 

X 

and of (i2{T;S), the second moment of the correlation function, 

^(T;S) = J2^(s(0)s(x)) c (6) 

X 

the "second-moment correlation length" £(T; S) is defined by 

The second field-derivative of the susceptibility Xa{T] S) is defined by 

X4 (T; S) = ^ (s(6)s(x)s(y)s(z)) c . (8) 

x,y,z 
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III. DEFINITIONS OF CRITICAL PARAMETERS 

In terms of the asymptotic behavior of these observables, we can now define the critical 
parameters, amplitudes and exponents that we are going to estimate using HT and LT series. 
The spontaneous magnetization has the asymptotic behavior 

M(")(T; S) ~ B(-\S)\t(S)f (s \l + a ( k \S)\t(S)\^ s) + ..) (9) 

as t(S) -> 0-. 

The asymptotic behaviour of the susceptibility as t(S) — > 0±, is expected to be 

X {±) (T; S) ~ C&\S)\t(S)\-TW(l + a£\S)\t(S)\ W + ... + b£\S)t(S) + . . .) (10) 
The correlation length 

£ {±) (T; S) ~ + 4^(5)1^)1^ + . . . + 6f ) (S , )t(S') + . . . ) (11) 

the specific heat 

C£\T; S)/k B ~ #'(5) In |f(S)|(l + ag^S)!^)^ + . . . + &g°(S)f(S) + . . .) (12) 
and the second field-derivative of the susceptibility x^[K\ S) 

X ( t\T; S) ~ -C| ±} (5) |t(S')|-^(l + at\S)\t(S)\ 6 ^ + ... + b?\S)t(S) + . . . ) (13) 

have analogous asymptotic behaviours. 

Different (universal) critical exponents (3(S), ^(S), v(S), ^(S) and different (non- 
universal) critical amplitudes B<-~\S), C^(S), f&(S) . . ., a^(S), af } (S), etc. are as- 
sociated with the various observables. We have reported in such detail our definitions of 
the critical amplitudes, because they differ significantly from those of other authors and it 
is necessary to use a consistent normalisation convention when comparing models expected 
to belong to the same universality class. Let us notice in particular that the estimates re- 
ported in the tables of Ref. 7 for the critical amplitudes of the susceptibility x^~\ u 'i are 
related to ours by the factor S 2 (-\nu c (S))' y / u c (S) . A similar remark applies to the specific 
heat amplitudes C^(u; S) for which the conversion factor is l/w c (S') 4 ' s '(lnu c (S')) 2 . Finally, 
the magnetization amplitudes of Ref. 7 should be multiplied by the factor (-ln-u c (S')) 1 / 8 / S to 
agree with ours. Of course, the amplitudes of the conformal field theory considered in the 
study of Ref. 17 are not comparable to our series quantities. 

As indicated in eqs.(9) - (13), for a given spin S, all asymptotic forms are moreover 
expected 18 to contain leading non-analytic confluent corrections characterized by the same 
exponent 9(S). Higher order corrections are also expected to contain logarithmic 18 factors. 
If universality holds, all exponents have to be S'-independent. 

The presence and the value of the confluent exponent has been discussed 19 ' 20 ' 21 ' 22,23 several 
times. From RG calculations 24 , both in the e-expansion and in the fixed- dimension approach, 
it was conjectured that 9 ~ 4/3 for the universality class of the two-dimensional Ising 
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model. Aharony and Fisher and later Blote and den Nijs argued 19 ' 22 that a^(S) = for 
S — 1/2 and indeed no such correction was revealed by the later very accurate study 13 ' 25 of 
the critical asymptotic expansion for x{K\ 1/2)- However, in the absence of more general 
results, the reliable assessment of the subleading asymptotic critical behaviour remained an 
open problem when S > 1/2. 



IV. ESTIMATES OF UNIVERSAL AMPLITUDE COMBINATIONS 

In terms of ^"^(If; S), ^ + \K;S) and X4~\k> S), a "hyper-universal" combination of 
critical amplitudes denoted by g( + \S) and usually called the "dimensionless renormalized 
coupling constant" , can be defined by 

Q (+)(S) = 3vC{ +) {S) 

Here the normalisation factor is chosen in order to match the usual field theoretic 

on 

definition 24 of g( + \S) and v denotes the volume per lattice site, measured in units of the 
square of a lattice constant. For all lattices one has v = a a 2 , with a the lattice constant. For 
the triangular lattice, a = v / 3/2, for the honeycomb lattice a = 3-^/3/4, and for the kagome 
lattice a = 2/^/3. 

We have also studied the hyper-universal combination usually denoted as 

R^\S) = (A^(S)/vy/ 2 (af^(S)) (15) 
and the Watson combination 26 

R C (S) = A( + \S)C (+ \S)/B(-\S) 2 (16) 
The other frequently considered universal combination 

R 4 (S) = Ci + \S)B^\S) 2 /C (+ \S) 3 , (17) 

is not independent of the previous ones, since Ra(S) = —^gl + \S)R^\S) 2 /Rc(S). 

All of these quantities are accurately known in the S = 1/2 case. As indicated in Ref. 14 ' 25 , 
it is known that #'(1/2) = |ln[tan(7r/8)] 2 « 0.4945385895, CW(l/2) ~ 0.962581732, 
and £(->(l/2) = 2 5 / 16 ln[l + v^] 178 ~ 1.22240995. In Refs. 23 ' 27 , we find the very accurate 
estimates Ci +) (l/2) = 4.379095(8) and / (+) (l/2) w 0.5670683. 

Therefore we can conclude # } (l/2) = 1.754364(2), #c(l/2) w 0.31856939, ^ +) (l/2) w 
0.39878194 and R 4 {l/2) = 7.336744(10). 

We have also considered the ratio C < - + \S)/C < -^' ) (S). In Ref. 25 , for S — 1/2, this ratio was 
computed with arbitrary precision to be C^(l/2)/C^(l/2) w 37.693652. 

Finally, we have estimated the ratio (S) / A^~\S) for various values of S. This ratio 
equals unity for S = 1/2 by self-duality. This was argued 28 in greater generality for the 
g-state (0 < q < 4) Potts model on the square lattice, which, for q = 2, reduces to the 
S — 1/2 Ising model. 

In what follows, we determine the values of these universal amplitude combinations for 
S > 1/2. The preliminary part of our series analysis is aimed at estimating the critical tem- 
peratures using the expansions of x^ + \K\ S) for S > 1. We employed a variety of methods: 
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Zinn- Justin improved-ratio formula 29 , Pade approximants (PA) and inhomogeneous differ- 
ential approximants (DA) 30 . The best results with DA's were obtained from approximants 
such that the polynomial coefficient of the highest derivative is even. (As a consequence, 
the approximants always contain an additional anti-ferromagnetic singularity at —K C (S), 
beside the one at K C (S)). Similarly to the LT analysis, but to a much smaller extent, the 
accuracy of our results tends to deteriorate with increasing S. In spite of this, our final 
HT estimates of the critical points, reported in Table III, show significant improvement in 
apparent accuracy and sizable discrepancies from the previous LT determinations 7 . 

For general values of S, less accurate estimates of K C (S) have been obtained in Ref. 31 from 
the ten term susceptibility series of Ref. 3 and from a renormalization group method. More 
recently other estimates 32 were obtained by a generalized cluster method. To our knowledge 
other accurate determinations of the critical points are available only for S = 1. They have 
been obtained either by analysing 21 the 21 term HT series of Ref. 6 for the susceptibility 
or by various transfer- matrix methods 33 ' 34 ' 35 . Some of these results have also been cited in 
Table III. 

We have then turned to the critical exponents j(S), v(S) and 74(5) and have evaluated 
them from the log-derivatives of the appropriate HT expansions by first order DA's biased 
with our HT estimates of the critical temperatures. This computation shows that the rel- 
ative variation of the exponents is smaller than 10 -3 , in the worst case, for S varying 
between 1/2 and 00. We report these results in Table IV without further details and simply 
conclude that universality and hyperscaling appear to be well supported for the leading 
critical exponents. 

It is perhaps also worth noticing that assuming the universality of 7 we can bias and 
therefore refine the determination of K C (S). This procedure does not change the central 
values of the critical points with respect to the unbiased one, but reduces the error bars. 

On the other hand, the estimate of the exponent 9(S) of the leading singular confluent 
corrections to scaling in the various observables remains quite elusive. Performing either a 
Baker-Hunter 36 or a Zinn- Justin 29 analysis, we can conclude that, at the level of accuracy 



TABLE III: Estimates of the critical inverse temperatures for the spin-S" Ising models on the 
sq lattice. Of course, the estimate Kc + \s), obtained from the HT series, must equal Kc \s) 
obtained from the LT series, and their common value is known exactly only for S = 1/2. For 
comparison, we have also reported other results beside those obtained from our HT and those 
obtained in Ref. 7 from the analysis of LT series. No error estimates are provided in Ref. 31 for the 
estimates of (S) obtained from the ten term series 3 of Camp and Van Dyke as well as for the 
estimates K C {S) obtained by a renormalization group method. 



S=l/2 


S=l 


S=3/2 


S=2 


S=5/2 


S=3 


S=oo 



K. 



0.44068679. . 
'~\S) 7 0.44068679.. 

0.441 



id +) (s) 31 

K C (S) 31 
K { C + \S) 33 
K C (S) 35 
K c {S) 3i 



0.458 



0.590473(5) 0.684255(6) 0.748562(8) 0.79541(1) 0.83106(2) 1.09315(2) 

0.5904727(9) 0.684338(46) 0.7487(14) 0.8025(35) 0.839(10) 

0.592 0.687 0.752 0.800 0.836 

0.610 0.704 0.770 0.818 0.855 

0.5904727(10) 
0.590471 
0.590076... 
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TABLE IV: Estimates of critical exponents obtained from our HT series for the spin-S" Ising 
models on the sq lattice. Of course the values of 7, v and 74, for S = 1/2 are exactly known. 





S=l/2 


S=l 


S=3/2 


S=2 


S=5/2 


S=3 


S=oo 




1.75 


1.7502(4) 


1.7500(4) 


1.7496(5) 


1.7500(5) 


1.7501(5) 


1.7494(8) 


u(S) 


1.0 


0.9999(6) 


0.9996(8) 


0.9994(8) 


0.9994(8) 


0.9994(8) 


0.9994(8) 




5.5 


5.498(4) 


5.497(5) 


5.497(5) 


5.497(5) 


5.497(5) 


5.497(5) 



made possible by the present extension of the HT series, the amplitudes of these corrections 
are very small, (or perhaps vanishing) for all values of S. We should mention that a similar 
conclusion was suggested for S — 1 in Ref. 22 , while the opposite conclusion was advocated 
in Ref. 21 . 

Once we have estimated the critical temperatures and verified the universality of the 
leading exponents, we can proceed with the analysis simply assuming that, for all values 
of S, these exponents take exactly the values expected for the universality class of the 
two-dimensional spin- 1/2 Ising model and using them along with our estimated critical 
temperatures to bias the evaluation of the HT and LT critical amplitudes defined by eqs.(9) 
- (13). 

Our estimates of the critical amplitudes are reported in Table V. 

We have employed quasi-diagonal non-defective PA's or DA's for extrapolating to K C (S) 
the effective amplitudes of the susceptibility and of the derivative of the specific heat, from 
the HT and the LT side of the critical points. We have similarly studied the effective ampli- 
tudes of the correlation-length (available only in the HT region) and of the magnetization. 
For proper comparison, in the same Table we have also cited the LT estimates of the crit- 
ical amplitudes for the spontaneous magnetization, the specific heat and the susceptibility 
previously obtained in Ref. 7 . These quantities have been multiplied by the above indicated 

TABLE V: Our series estimates of the critical amplitudes on the HT and the LT side of the critical 
point for various spin-S" Ising models on the square lattice. Along with our estimates of the LT 
amplitudes we have reported for comparison also the results of Ref. 7 multiplied by the proper 
conversion factors. The values of C ( -+\l/2), A ( - + \l/2), C { ~\\/2), BH(l/2) and #)(l/2) are 
exactly known, those of /( + )(l/2) and c{ + \l/2) are very accurately known. 





S=l/2(ex.) 


S=l 


S=3/2 


S=2 


S=5/2 


S=3 


S=oo 


C<-+\S) 


0.9625817... 


0.5514(2) 


0.4307(2) 


0.3755(2) 


0.3441(2) 


0.3254(3) 


0.2351(2) 


A(+\S) 


0.494538589 . . . 


0.736(10) 


0.854(4) 


0.917(4) 


0.956(4) 


0.983(4) 


1.054(6) 


/(+)(£) 


0.567068. . . 


0.4640(2) 


0.4309(2) 


0.4159(2) 


0.4082(2) 


0.4033(2) 


0.3900(2) 


ci + \s) 


4.379095(8) 


0.9630(4) 


0.5073(3) 


0.3591(3) 


0.2902(3) 


0.2533(3) 


0.1239(3) 


C7(-)(5) 


0.02553697... 


0.01462(3) 


0.0114(3) 


0.0102(6) 


0.0090(8) 


0.0055(30) 








0.01462(2) 


0.0109(29) 


0.0094(10) 


0.0096(33) 








1.22240995... 


1.131(4) 


1.076(5) 


1.041(5) 


1.016(5) 


1.001(5) 




b(-\s) 7 




1.1313(2) 


1.077(9) 


1.042(16) 


1.030(19) 


1.016(26) 




A(~\S) 


0.494538589 . . . 


0.738(6) 


0.855(10) 


0.915(10) 


1.1(2) 


1.1(2) 




A(-\S) 7 




0.73(2) 


0.77(6) 


0.86(8) 


0.87(9) 
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TABLE VI: Universal combinations of critical amplitudes for various spin-5 Ising models on the 
square lattice. The exactly (or very accurately) known values for S = 1/2 are reported in the first 
column. For S > 1/2, the series estimates of this note are reported in the successive columns. In 
the last line we have reported the estimates of Rc(S) obtained by combining our present estimates 



of #'(5) and C<- 


^~\S) with the estimates of \S) g 


iven in Ref. 7 . 








S=l/2(ex.) 


S=l 


S=3/2 


S=2 


S=5/2 


S=3 


S = oo 


c(+\s)/c(-\s) 


37.693652 . . . 


37.71(9) 


38(1). 


37(2). 


38(3). 


59(32). 




A(+\S)/A(-\S) 


1.0 


0.997(21) 


0.999(16) 


1.0(1) 


0.87(16) 


0.89(17) 




9 P\S) 


1.754364(2) 


1.753(2) 


1.753(2) 


1.752(3) 


1.752(3) 


1.752(3) 


1.752(3) 


R ( +\S) 


0.39878194... 


0.398(3) 


0.398(1) 


0.398(1) 


0.399(1) 


0.400(2) 


0.400(2) 


Rc(S) 


0.31856939... 


0.317(5) 


0.318(2) 


0.318(2) 


0.319(3) 


0.319(5) 




Rc(S) 7 




0.317(4) 


0.317(7) 


0.317(11) 


0.31(1) 


0.31(1) 





conversion factors to agree with our normalisation conventions. The uncertainties we have 
reported, which allow for the observed spreads in the approximant values, provide a sub- 
jective assessment of residual trends in the sequence of estimates and for the (unbiased) 
uncertainties of the critical points. The HT amplitudes can be determined with a relative 
accuracy ranging from « 10~ 3 in the case of the susceptibility, to ps 1CT 2 in the case of the 
specific heat. The LT amplitudes are subject to larger relative uncertainties, increasing with 
S, and reaching up to tv 50% for S > 2. In some cases, in order to improve the accuracy 
of the estimates of the LT amplitudes for S < 2, we have based our extrapolations only on 
the data for \t(S)\ ^ 0.02 — 0.04. This unconventional but reasonable procedure reduces the 
sensitivity of the approximants to the unphysical nearby singularities. Unfortunately, even 
this prescription fails to work satisfactorily for S > 2. 

Using only the HT series, we can evaluate g( + \S), either directly, in terms of the am- 
plitudes reported in Table V, or by extrapolating to the critical points via DA's the HT 
expansion of the inverse effective coupling l/g( + \K; S) = —^-^(K; S) 2 x(K; S) 2 /xa{K\ S). 
A third approach consists in studying the residua at x — 1 of the series with coefficients 
a n (S) = c n (S)/d n (S), where c n (S) are the HT coefficients of £ 2 (K;S) and d n (S) are the 
coefficients of the quantity Xi{K\ S)/x 2 (K'i S). In Table V we have reported the results of 
the latter procedure since it yields estimates with smaller spreads. 

Several other estimates 23 ' 37 ' 38 ' 39 ' 40 ' 41 obtained by a variety of methods are also available 
in the literature. 

Using also the LT series, we have evaluated, directly in terms of the single amplitudes, 
the other mentioned universal combinations, for a range of values of S. We have reported in 
Table VI, our series estimates of all these quantities for S > 1/2. In conclusion, whenever 
only HT amplitudes are involved, our estimates, within a precision up to 0.1%, are indepen- 
dent of S, in full agreement with universality. On the other hand, our reanalysis of the LT 
series has been only partially successful: whenever LT amplitudes also enter into the com- 
binations, universality appears to be fairly well respected for S < 5/2, but the uncertainties 
grow notably larger for larger values of the spin. 

In the next section we describe the theory of lattice-lattice scaling, and show how it can 
be used to extend our estimates of the critical amplitudes from the square lattice to other 
two-dimensional lattices. 
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V. LATTICE-LATTICE SCALING 



The theory of lattice-lattice scaling was developed by Betts, Guttmann and Joyce 15 in 
the early '70s. It explains how amplitudes change within a given universality class, as one 
moves from one lattice to another. It can also be viewed as a generalisation of the law of 
corresponding states. In this section we give a terse development of the theory, and apply 
it to the problem at hand. 

In order to review the general ideas let us first consider the Weiss theory or mean field 
theory of a magnetic system. The equation of state is well known to be 

h= -m 3 {l + 3t/m 2 ). 
3 

Here t = T/T c - 1, h = fiH/kT, and m = M(T)/M(0) are the reduced temperature, 
magnetic field and magnetization, respectively. 

Then the law of corresponding states says that the equation of state is the same for all 
lattices. That is, 

mx(t, h) = my(f, h), 

where X and Y denote two lattices. That is to say, the lattice dependence is entirely 
contained in the critical temperature T c . 

A more complex model is the three-dimensional spherical model, for which the critical 
equation of state is: 

h = D x m 5 {l + t/m 2 ) 2 . 
Here both T c and the amplitude D are lattice dependent. Thus 

m x (t, hx) = m Y (t,h Y ). 

We see that we must scale the field variable, so that jfc = but that there is no need to 
scale the reduced temperature. 

Let us now consider the case of the (zero- field) spin S — 1/2 Ising model on the triangular 
(T) and hexagonal (H) lattices. The star-triangle relation 42 ' 43 allows us to relate the free- 
energy, susceptibility and spontaneous magnetization between the lattices: 

2f H (K H ) = h(K T ), 
M Ti0 (*t) = M Hl o(t H ), 
2xt(v t ) = Xh(^h) + Xh(-uh), 

where / = 1/NlnZ, Kx = Jx/kT, and fx = tanh(Jx/A;T). Here we see that the reduced 
temperature needs to be re-scaled for the free-energy to be universal. This is not restricted 
to the triangular-honeycomb pair, but in that case it is easy to be explicit. 

All these examples can be encapsulated in the following expression for the singular part 
of the free-energy: 

nxfx(tx, h x ) = n Y f Y (tY, h Y ) = f(t, h), 
where the reduced temperature and field are scaled by 

nxhx = "^y^y = h 

and 

gxtx = gyt Y = t. 
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The singular part of the free-energy, f(t, h), is then a universal (lattice independent) function 
for a given model. 

Equivalently, by differentiation we obtain 

m x (tx, hx) = m Y (t Y , h Y ) = m(t, h), 

and 

Xx(t x , hx)/nx = Xv(t Y , h Y )/n Y = x(t, h), 
where m(t, h) and \{t, h) are universal functions for a given model. 
Writing mx = B^\—txY ■, it follows that 

9x = I Ux_ 

9y \ yJ B Y " ) 

Using this result, and the exact scaling parameters gx and g Y given below, it is a trivial 
matter to calculate the magnetization amplitudes for the other lattices we consider (trian- 
gular, hexagonal and kagome (K)), taking as input the square lattice amplitudes given in 
Table V. These amplitude estimates are given in Tables VII, VIII and IX. 
Writing xx = Qx^x 7 ^ follows that 

4+74+) = i* (^V 7 = 4-)/4-). 

n Y \g Y J 

Similarly, it is a trivial matter to calculate the susceptibility amplitudes for the other lattices 
we consider, taking as input the S — 1/2 square lattice amplitudes given in Table V. 

Further differentiation gives the corresponding relationship for higher field derivatives, 
and we readily obtain 

for the high-temperature field derivatives, (where only the even-order derivatives are non- 
zero). The corresponding result for low-temperature field derivatives is 

Taking temperature derivatives, one readily establishes that the specific heat amplitudes 
satisfy 

4 +, /4 +, = ^(^y~* = 4- ) M- ) . 

nx \9y j 

As a = for the two-dimensional Ising model, this simplifies to 

4 + >/4+> = ^(»y = 4->/4->. 

n x \9y j 

We have similarly calculated the specific-heat amplitudes for the other lattices we consider, 
taking as input the square lattice amplitudes given in Table VI. These are also given in 
Tables VII, VIII, and IX. 
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TABLE VII: Estimates of the critical amplitudes on the HT and the LT side of the critical point 
for various spin-5 Ising models on the triangular lattice, as obtained by lattice-lattice scaling, 
using the square- lattice series estimates as a basis. The values of CW(l/2), A ( - + \l/2), C^\l/2), 
£(-)(l/2) and iH(i/2) are exactly known. 





S=l/2(ex.) 


S=l 


S=3/2 


S=2 


S=5/2 


S=3 


S=oo 


C(+)(5) 


0.92420696 . . . 


0.5294(2) 


0.4135(2) 


0.3605(2) 


0.3304(2) 


0.3124(3) 


0.2257(2) 


A(+\S) 


0.499069377 . . . 


0.743(10) 


0.862(4) 


0.925(4) 


0.965(4) 


0.992(4) 


1.064(6) 




0.525315... 


0.4298(2) 


0.3992(2) 


0.3853(2) 


0.3781(2) 


0.3736(2) 


0.3613(2) 


cl +) (5) 


4.000248(8) 


0.8797(4) 


0.4634(3) 


0.3280(3) 


0.2651(3) 


0.2314(3) 


0.1132(3) 




0.024518902 . . . 


0.01404(3) 


0.0109(3) 


0.0098(6) 


0.0086(8) 


0.0053(30) 






1.203269903... 


1.113(4) 


1.059(5) 


1.025(5) 


1.000(5) 


0.985(5) 






0.4990693724 . . . 


0.745(6) 


0.863(10) 


0.923(10) 


1.1(2) 


1.1(2) 





TABLE VIII: Estimates of the critical amplitudes on the HT and the LT side of the critical point 
for various spin-S* Ising models on the hexagonal lattice, as obtained by lattice-lattice scaling, 
using the square-lattice series estimates as a basis. The values of CW(l/2), A^(l/2), CH(l/2), 
S( - )(i/2) and iH(i/2) are exactly known. 





S=l/2(ex.) 


S=l 


S=3/2 


S=2 


S=5/2 


S=3 


S=oo 




1.0464170... 


0.5994(2) 


0.4682(2) 


0.4082(2) 


0.3741(2) 


0.3537(3) 


0.2556(2) 


A(+\S) 


0.4781063817... 


0.712(10) 


0.826(4) 


0.887(4) 


0.924(4) 


0.950(4) 


1.019(6) 


/(+)(£) 


0.657331. . . 


0.5379(2) 


0.4995(2) 


0.4821(2) 


0.4732(2) 


0.4675(2) 


0.4521(2) 


C[ +) (S) 


5.352965(8) 


1.1772(4) 


0.6201(3) 


0.4390(3) 


0.3547(3) 


0.3096(3) 


0.1515(3) 


c ( -\s) 


0.027761095 . . . 


0.01589(3) 


0.0124(3) 


0.0111(6) 


0.0098(8) 


0.0060(30) 




B(~\S) 


1.253177691... 


1.159(4) 


1.103(5) 


1.067(5) 


1.042(5) 


1.026(5) 






0.4781063817... 


0.714(6) 


0.827(10) 


0.885(10) 


1.1(2) 


1.1(2) 





TABLE IX: Estimates of the critical amplitudes on the HT and the LT side of the critical point for 
various spin-5 Ising models on the kagome lattice, as obtained by lattice-lattice scaling, using the 
square-lattice series estimates as a basis. The values of C(+)(l/2), #'(1/2), C ( -~\l/2), B( _ )(l/2) 
and A^~\l/2) are exactly known. 





S=l/2(ex.) 


S=l 


S=3/2 


S=2 


S=5/2 


S=3 


S=oo 


CW(S) 


1.01814223... 


0.5832(2) 


0.4556(2) 


0.3972(2) 


0.3640(2) 


0.3442(3) 


0.2487(2) 


A(+\S) 


0.4800615653 . . . 


0.714(10) 


0.829(4) 


0.890(4) 


0.928(4) 


0.954(4) 


1.023(6) 


/«(£) 


0.618474. . . 


0.5061(2) 


0.4700(2) 


0.4536(2) 


0.4452(2) 


0.4399(2) 


0.4254(2) 


c[ + \s) 


5.046953(8) 


1.1099(4) 


0.5847(3) 


0.4139(3) 


0.3345(3) 


0.2919(3) 


0.1428(3) 


c<--\s) 


0.02701097... 


0.01546(3) 


0.0121(3) 


0.0108(6) 


0.0095(8) 


0.0058(30) 




B(~\S) 


1.238655888... 


1.146(4) 


1.090(5) 


1.055(5) 


1.030(5) 


1.014(5) 




A(~\S) 


0.4800615653 . . . 


0.716(6) 


0.830(10) 


0.888(10) 


1.1(2) 


1.1(2) 
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Finally, the correlation function amplitudes were calculated exactly from various star- 
triangle transformations by Thompson and Guttmann 44 for the true correlation length. It 
follows from universality that these same transformations should hold also for the second- 
moment correlation length amplitudes. We show below that this is equivalent to the uni- 
versality of R^~\S), which is a conclusion of this work. 

In 44 it was shown by explicit calculation that 

/W(l/2)^(l/2)//«(l/2)^(l/2) = 1/3. 

# ) (1/2)< + s ) (1/2)//W(1/2)kW(1/2) = Vs. 
/i +) (l/2)K«(l/2)//4 +) (l/2)< + H ) (l/2) = 2/3. 

The universality of R^~\S) = ( ) a f + (S), taken together with the above lattice- 
lattice scaling relation for the specific heat amplitude A + (S) implies the lattice-lattice scaling 
relation 

A+)/ ,w = r 

9x \cr Y n Y J 

where a is the area per site, and ax = 1, V / 3/2, 3a/3/4 and 2/y/3 for the square, triangular, 
honeycomb and kagome lattices respectively. This is equivalent to the explicit amplitude 
scaling reported in 44 , and confirms the expectation that the true correlation length amplitude 
and the second-moment correlation length amplitude scale similarly. 

For the 2d Ising model we can calculate the scaling parameters n x and gx exactly from 
known spin-1/2 spontaneous magnetization and specific heat amplitudes. The critical points 
are also exactly known. These are given in the table below: 

X K c n x gxK c ,T/K c ,x 

Tr. ln(3)/4 1 1 

Sq. -m(y2-l)/2 3^3/4 1/V2 

Ho. -\n(2-V3)/2 2 l/y/3 

Ka. -]n(2/«s/3-l)/4 S^v^-l) 16 9( v / 3-l) 8 



In addition to the above results, one can also derive scaling relations for the amplitudes 
of sub-dominant singularities. For example, consider the susceptibility of the 2d Ising model 
on lattice X. Writing 

Xx ~ C*t-~< + Cft 1 ^, 
from lattice-lattice scaling we can derive the following amplitude relations: 

9sL = nx ( 9Y 
Co n Y \gx 

Cl = nx(9YV-\ 
Ci n Y \gx) 

The second expression is false for the kagome lattice. It is corrected by the theory of extended 
lattice-lattice scaling developed by Gaunt and Guttmann 16 in 1978. In that theory, a third 
scaling parameter needs to be introduced. 
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It is widely accepted, and in complete agreement with the results of the first part of this 
paper, that the spin-5 1 Ising model is in the same universality class as the spin-1/2 Ising 
model. This is the only assumption we require in order to apply the theory of lattice-lattice 
scaling to the square lattice data given in the previous section. The only subtlety is whether 
"universality" really extends to lattice-lattice universality. While this assumption seems 
natural, we did attempt to verify it by estimating amplitudes for other lattices from the 
rather short data available in 3 . The longest effective series is the triangular lattice series. We 
found the high-temperature susceptibility amplitude, as estimated by Pade approximants, 
to be C( + )(l) = 0.529 from this series, in complete agreement with the more precise value 
C (+) (l) = 0.5294(2) found by lattice-lattice scaling, and reported in Table VII. 

Note too that there is no loss of accuracy, as the conversions from lattice to lattice are 
exact. For example, based on the recent estimate of the leading susceptibility amplitude of 
the spin 1/2 square lattice Ising model given in 13 , application of lattice-lattice scaling gives 
the corresponding amplitude on the triangular lattice to the same precision, viz: 

Cl = 0.9242069582451643296971575778559317176696261520028389. 

Similarly accurate results for other lattices can be readily written down, as can equally 
accurate sub-dominant amplitudes. 
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